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With the superconducting cuprates in mind, a set of unitary transformations was used to decouple 
electrons and phonons in the strong-coupling limit. While phonons remain almost unrenormalised, 
electrons are transformed into itinerent singlet and triplet bipolarons and thermally excited polarons. 
The triplet/singlet exchange energy and the binding energy of the bipolarons are thought to account 
for the spin and charge pseudogaps in the cuprates, respectively. We calculated the Hall Lorenz 
number of the system to show that the Wiedemann- Franz law breaks down due to the interference 
of the polaron and bipolaron contributions to heat flow. The model provides a quantitative fit to 
magnetotransport data in the cuprates. Furthermore we are able to extract the phonon component 
of the thermal conductivity with the use of experimental data and the model. Our results further 
validate the use of a charged Bose gas model to describe normal and superconducting properties of 
unconventional superconductors. 



I. INTRODUCTION 

The theory of high temperature superconductivity 
is one of the biggest challenges in condensed mat- 
ter physics today and is reflected by the large num- 
ber of the theoretical modelsi* 2 * 3 ^*^ proposed to date. 
Most of these theories have a foundation built on the 
unconventional Cooper pairs whilst others concentrate 
on the basic principals of the symmetry breaking and 
competing orders. On the other hand the electron- 
phonon interaction continue to gather support through 
isotope effect measurements^, infrarec&iLiii and thermal 
conductivity!!, neutron scattering^, and more recently 
in ARPESi&il To account for the high values of T c 
in the cuprates, it is necessary to have electron-phonon 
interactions larger than those found in the intermedi- 
ate coupling theory of superconductivity 1 ^. Regard- 
less of the adiabatic ratio, the Migdal-Eliashberg theory 
of superconductivity and Fermi-liquids has been shown 
to breakdown at A = 1— using the (1/A) expansion 
technique^. The many-electron system collapses into the 
small (bi)polaron regimaiSiiSiiS at A > 1 with well sepa- 
rated vibration and charge-carrier degrees of freedom. 

Remarkably the discovery of high T c cuprates was mo- 
tivated by the polaron model22i. The isotope effect, high 
values of the static dielectric constants and optical spec- 
troscopy in the cuprates suggest that the electron-phonon 
interaction is sufficient to bind small polarons into small 
bipolarons. At first sight these carriers have a mass too 
large to be mobile, however it has been shown that the 
inclusion of the on-site Coulomb repulsion leads to the 
favoured binding of intersite oxygen holes 2,21 . The inter- 
site bipolarons can then tunnel with an effective mass of 
only 10 electron masse o 21 i 22 ' 23 i 24 . 

Mott and one of the authors (ASA) proposed a simple 
model^ of the cuprates based on bipolarons, Fig.l. In 
this model all the holes (polarons) are bound into small 
intersite singlet and triplet bipolarons at any tempera- 



ture. Above T c this Bose gas is non-degenerate and below 
T c phase coherence (ODLRO) of the preformed bosons 
sets in, followed by superfluidity of the charged carriers. 
The triplet and singlet states are separated by an ex- 
change energy J which explains the spin gap observed in 
many NMR and neutron scattering experiments^^!. Of 
course, there are also thermally excited single polarons 
in the model. Their density becomes comparable with 
the bipolaron density at the temperature T* which is 
about half of the bipolaron binding energy A. There 
is much evidence for the crossover regime at T* and 
normal state pseudogap in the cuprates 2 ^. For exam- 
ple, measurements of single particle tunneling and An- 
dreev reflections 2 - clearly demonstrate that there exists 
two distinct energy gaps. These energy gaps are half 
of the binding energy of preformed bipolarons and the 
coherence energy gap of the superconducting phase as 
proposed in Refill. 

Other experimental observations which have been 
satisfactorily explained using this particular approach 
include Hall ratio^!, Hall angle^i, in-plane resistivity^, 
infrared conducitivity^, magnetic susceptibility**, 
tunneling spectroscopy 3 **, isotope effect, thermal 
conductivity*^ and the upper critical field*! in the 
superconducting state. ARPES measurements indicate 
the presence of an angular dependent narrow peak and 
a featureless background. In the polaronic model the 
ARPES spectrum can be explained if one considers a 
charge transfer Mott-insulator and the single polaron 
spectral function^. The conclusion is that the ARPES 
spectrum can in fact be explained with the absence or 
with a small Fermi surface. Specific heat measurements 
have a striking resemblance to the superfluid transi- 
tion of weakly-interacting 3D bosons 3 -. A theoretical 
models, experimental data and analysis of the specific 
heat in a magnetic field*! suggest this is a convincing 
argument. 

In this paper we develop the theory of the normal state 
magnetotransport in bipolaronic systems and explain the 
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FIG. 1: Bipolaron picture of high temperature superconduc- 
tors. A corresponds to the singlet intersite bipolaron. B is 
the triplet intersite bipolaron, which naturally includes the 
addition of an extra excitation band. The crosses are copper 
sites and the circles are oxygen sites, w is the half bandwidth 
of the polaron band, t is the half bandwidth of the bipolaron 
band, A/2 is the bipolaron binding energy per polaron and J 
is the exchange energy per bipolaron. 



netic field. As a result, the "Hall" Lorenz number, 
Lh = (e/kB) 2 K xy / l {Ta xy ) 1 has been directly measured 
in Y Ba2CusOe.g5 because transverse thermal n xy and 
electrical a xy conductivities involve only the electrons. 
The experimental L xy showed a linear temperature de- 
pendence, which violates the Wiedemann-Franz law. It 
is clear that it would be difficult to explain this experi- 
mental observation in the framework of any Fermi-liquid 
model. 

A charge 2e Bose gas of bipolarons naturally leads to 
the conclusion^ that the Lorenz number in the normal 
state is strongly suppressed due to the double charge 
per carrier and classical statistics of bipolarons com- 
pared with the Lorenz number in the normal metal, 
Lq = 7r 2 /3. Here we demonstrate that the Wiedemann- 
Franz law breaks down because of the interference of the 
polaron and bipolaron contributions to the heat trans- 
port. When thermally excited polarons and triplet pairs 
are included, the bipolaron model explains the violation 
of the Wiedemann-Franz law in the cuprates and the Hall 
Lorenz number as seen in the experiments. 



II. LOW FERMI ENERGY: PAIRING IS 
INDIVIDUAL IN MANY CUPRATES 



experimental Hall Lorenz number for Y BaiCu^O^+x as 
measured by Zhang et al4^. The case of degenerate sin- 
glet and triplet bipolarons was briefly discussed by us in 
Refi^. Particular interest within this paper lies in the 
conclusion that the Wiedemann-Franz law is violated in 
cuprates. This departure from the Fermi liquid picture is 
seen in both the superconducting and normal state and 
might be related by a common mechanism 4 ^. 

The extraction of the electronic thermal conductivity 
term has proven difficult as both the electronic term, K e i 
and the phononic term, K p h make a comparable contri- 
bution. However recent experiments appear to have got 
around this problemS 4 ^^. Takenaka et al4^ found that 
the insulating state thermal conductivity (therefore pre- 
dominantly phononic) was approximately independent 
of doping except at y ~ 1 — 0.75 in Y BaiCu-zOi-y. 
Therefore the phononic term for the thermal conductiv- 
ity above y ~ 0.5 could be estimated. This analysis led 
to the conclusion that the K e i is constant or weakly T- 
dependent in the normal state like that found by Salamon 
et al^. This approximately T-independent K e i in the 
underdoped region therefore implies the violation of the 
Wiedemann-Franz law (since the resistivity is found to 
be a non- linear function of temperature in this regime). 
This breakdown of the Wiedemann-Franz law has been 
seen in other cuprates such as Pri- x Ce x CuO^ at op- 
timal doping 4 ^ and Lai- x Sr x CuO^ at undcrdoping 47 . 
Zhang et al* recently developed a novel method to de- 
termine the Lorenz number, based on the thermal Hall 
conductivity. The thermal Hall effect is a purely elec- 
tronic effect since phonons are not affected by the mag- 



The possibility of real-space pairing, as opposed to 
BCS-like pairing, has been the subject of much discus- 
sion. Experimental and theoretical evidence for an ex- 
ceptionally strong electron-phonon interaction in high 
temperature superconductors is now so overwhelming, 
that even advocates of nonphononic mechanisms 4 ^ ac- 
cept this fact. Nevertheless the same authoraS^ dis- 
miss any real-space pairing claiming that pairing is col- 
lective in cuprates. They believe in a large Fermi surface 
with the number of holes (1 + x) rather than x in super- 
conducting cuprates, where x is the doping level like in 
La2- x Sr x CuOi. As an alternative to a three-dimensional 
Bose-Einstein condensation of bipolarons these authors 
suggest a collective pairing (i.e the Cooper pairs in the 
momentum space) at some temperature T* > T c , but 
without phase ordering. In this concept the phase co- 
herence and superconducting critical temperature T c are 
determined by the superfluid density, which is propor- 
tional to doping x due to a low dimensionality, rather 
then to the total density of carriers (1 + x). On the 
experimental side a large Fermi surface is clearly incom- 
patible with a great number of thermodynamic, mag- 
netic, and kinetic measurements, which show that only 
holes doped into a parent insulator are carriers in the 
normal state. On the theoretical ground this preformed 
Cooper-pair (or phase-fluctuation) scenario contradicts a 
theorem 4 ^, which proves that the number of supercarri- 
ers (at T = 0) and normal-state carriers is the same in 
any clean superfluid. 

Objections against real-space pairing also contradict a 
parameter-free estimate of the Fermi energy, which shows 
that the pairing is individual in cuprates^. Indeed, first- 
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principles band structure calculations show that copper, 
alkali metal, and magnesium donate their outer electrons 
to oxygen, Cqo, and boron in cuprates, doped fullerenes, 
and in M<?i?2, respectively. In cuprates the band struc- 
ture is quasi-two-dimensional (2D) with a few degenerate 
hole pockets. Applying the parabolic approximation for 
the band dispersion we obtain the renormalized Fermi 
energy as 



(1) 



where d is the interplane distance, and n,, m* arc the den- 
sity of holes and their effective mass in each of the hole 
subbands i renormalized by the electron-phonon (and 
by any other) interaction. One can express the renor- 
malized band-structure parameters through the in-plane 
magnetic-field penetration depth at T ~ 0, measured ex- 
perimentally: 
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As a result, we obtain a parameter-free expression for the 
"true" Fermi energy as 



h 2 dc 2 



(3) 



where g is the degeneracy of the spectrum which might 
depend on doping in cuprates. One expects 4 hole pock- 
ets inside the Brillouin zone (BZ) due to the Mott- 
Hubbard gap in underdoped cuprates. If the hole band 
minima are shifted with doping to BZ boundaries, all 
their wave vectors would belong to the stars with two 
or more prongs. The groups of wave vectors for these 
stars have only ID representations. It means that the 
spectrum will be degenerate with respect to the number 
of prongs which the star has, i.e g ^ 2. Because Eq.(3) 
does not contain any other band-structure parameters, 
the estimate of cf using this equation does not depend 
very much on the parabolic approximation for the band 
dispersion. 

Generally, the ratios n/m* in Eq.(l) and Eq.(2) are not 
necessarily the same. The 'superfluid' density in Eq.(2) 
might be different from the total density of delocalized 
carriers in Eq.(l). However, in a translationally invari- 
ant system they must be the samA This is also true 
even in the extreme case of a pure two-dimensional su- 
perfluid, where quantum fluctuations are important. One 
can, however, obtain a reduced value of the zero tempera- 
ture superfluid density in the dirty limit, I <C £(0), where 
£(0) is the zero-temperature coherence length. The latter 
was measured directly in cuprates as the size of the vor- 
tex core. It is about 10 A or even less. On the contrary, 
the mean free path was found surprisingly large at low 
temperatures, I ~ 100-1000 A. Hence, it is rather prob- 
able that novel superconductors are in the clean limit, 
I £(0), so that the parameter- free expression for ep, 
Eq.(3) holds. 



TABLE I: The Fermi energy (multiplied by the degeneracy) 
of cuprates 

Compound T c (K) \ H , ab (A) d(A) ge F (meV) 
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Parameter-free estimates of the Fermi energy obtained 
by using Eq.(3) are presented in the Table. The renor- 
malised Fermi energy in more than 30 cuprates is defi- 
nitely less than 100 meV. 

In many cases (Table) the renormalized Fermi energy is 
so small that pairing is certainly individual. Such pairing 
will occur when the size of a pair, p is smaller than the 
inter-pair separation, r. The size of a pair is generally 



P : 



(4) 



where A is its binding energy. The separation of pairs 
can be directly related to the Fermi energy in 2D 



= h 



(5) 
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We see that the true condition for real-space pairing is 

£f<7tA. (6) 

The bipolaron binding energy is thought to be twice 
the so-called pseudogap-. Experimentally measured 
pseudogap of many cuprates^S is as large as A/2 > 
50me V, so that Eq. (6) is well satisfied in underdoped and 
even in a few optimally and overdoped cuprates. One 
should notice that the coherence length in the charged 
Bose gas has nothing to do with the size of the boson as 
erroneously assumed by some authors^. It depends on 
the interparticle distance and the mean-free path, and 
might be as large as in the BCS superconductor. Hence, 
it would be incorrect to apply the ratio of the coherence 
length to the inter-carrier distance as a criterium of the 
BCS-Bose liquid crossover. The correct criterium is given 
byEq.(6). 



where we have set h = 1 and also from now on set ks — 
c = 1. v = dE/dli, t is the relaxation time and we as- 
sume that it depends on the kinetic energy, E = k 2 /2m. 
F = (E - n) S7T/T + V(/j - 2e0) and g = g s = 2e/m s for 
singlet bipolarons with the energy E = k 2 /{2m s ). For 
triplet bipolarons, F = (E + J - ^)VT/T + V(/x - 2e0), 
g = g t = 2e/m t and the energy E = k 2 /(2m t ). F = 
(E + A/2 - n/2)VT/T + V(/x/2 - e<f>), E = k 2 /{2m p ) 
and g = g p — e/m p for thermally excited polarons. Here 
m Sj t. p are the singlet and triplet bipolaron and polaron 
masses of two-dimensional carriers. 

Eqns.(7-ll) are used to find the thermal and electrical 
currents induced by the applied thermal and potential 
gradients in a magnetic field. 



IV. ELECTRICAL AND THERMAL CURRENTS 



III. NON-EQUILIBRIUM BIPOLARON AND 
POLARON DISTRIBUTION FUNCTIONS 

Because thermally excited phonons and (bi)polarons 
are well decoupled in the strong-coupling regime of 
the electron-phonon interaction the standard Boltzmann 
equation for kinetics of renormalised carriers may be ap- 
plied. The distribution function for each carrier is given 
as 



/ = fo + h 



(J) 



Here fo is the distribution function at equilibrium and 
fi (k) is the deviation of the distribution function away 
from equilibrium. We assume that fx is small compared 
t° fo- fo (polaron) is the Fermi-Dirac distribution and 
both /q and fo (singlet and triplet bipolaron, respec- 
tively) are the Bosc-Einstein distributions, 
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(8) 



(9) 
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fi is the chemical potential, J is the exchange energy 
which separates the triplet state from the singlet state, 
and A is the bipolaron binding energy per pair which is 
assumed to be of s-symmetry. The latter assumption of 
s-wave bulk pairing symmetry of a single-particle gap has 
been shown to be a valid one22^iiS. 

We make use of the r— approximation^ in the pres- 
ence of the electric field E, temperature gradient VT and 
magnetic field B|| z _L E and VT, 



/(k) = / (£)+T-^ 
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The electrical current for each carrier is given by 

r=g£v/f(k) (12) 

k 

where a = s,p, t, and q is the carrier charge. We find the 
x direction component as 

j° = a^V a (/z-2e^)+< y V y (M-2e«£) 

+b% x V x T + b% y V y T (13) 

and also the y direction component 
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and 



where 



, r . = / °° dEEr r jE)[\ + (g a r a (E)Bf]^dfS/dE 
{Ta) S^dEfZ 

The number densities, n a of the three carriers can be 
evaluated as 
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and the results are shown in Fig. 2 for optimally doped 
Y BCL2CU3O6.95 using experimental estimates of the en- 
ergy gaps as measured by Mihailovic et al2&. 




FIG. 2: number density ratios for triplets and polarons in the 
optimally doped regime 

To find the thermal currents we calculate the energy 
flux for each carrier which are given by 



(20) 



where EP = E + A/2 + e<j>, E s = E + 2ecj) and E t 
E + J + 2e<f>. The x direction component is 
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(23) 



V. TRANSPORT COEFFICIENTS 

In thermal transport experiments one has to make sure 
the electrical current is set to zero j = j P + j s + jt = 0. 
Using eqns. (13) and (14), we find, 
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where b xx = b p x + b xx + b xx , a yx = a p x + af,„ + al.„, etc. 



By substituting eqn (24) and (25) into eqn (21) and (22) 
we obtain for the thermal conductivities k xx , and for the 
thermal Hall conductivity, n yx 
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and 
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The electrical conductivity is denned in the presence 
of an electric field E = — V</> and in the absence of a 
thermal and chemical gradient, VT — V/i = 0, 
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We can combine cqns. (26), (27), (30) and (31) to obtain 
the Lorenz number and the Hall Lorcnz number as 
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respectively. 

We also define the Hall ratio. The Hall effect is the 
appearance of an electric field perpendicular to both the 
current and applied magnetic field, which is perpendicu- 
lar to the current. If the current flows in the x-direction 
only, j x = j and j y = 0, then we can eliminate E x and 
find the Hall ratio, which is defined as Rh = E y /(Bj), 
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VI. WEAK-FIELD APPROXIMATION 

If the magnetic field is weak, g a T a B « 1, we can 
ignore all terms in B 2 and higher order. The previous 
definition of (r£) simplifies and becomes 
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Next we can apply the coefficients, eqns (15) and (23), 
into the transport coefficients derived in the previous sec- 
tion. 



The resistivity, Hall resistivity and Hall ratio for a 
triplet, singlet bipolaron and polaron system are respec- 
tively 
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A 



s2,t2 



_ n Si t(T Sit }m p 
n p (T p )m s , t 

9s,tn s ,t(T^ t )m p 
9 P n P (T 2 }m Sv t 



(36) 
(37) 



(38) 



(39) 



(40) 



With the use of eqns. (32) and (33) we arrive at the 
Lorenz number 

L = D 1 [L p + AA sl L s + 4A n L t ] + (41) 

£?L4 sl (2r p -r s + A/T) 2 + 
A tl (2T p -T t + (A - J)/T) 2 + 

4A sl A tl (r t ~r s + j/T) 2 }} 

and the Hall Lorenz number 

L H = D 2 (L P + 4A s2 L s + AA t2 L t ) + D X D 2 x (42) 
[(4A 2 sl + 2A s2 )(A/T + 2T P - T s ) 2 + 
4A 2 tl ((J 2 + A 2 )/T 2 + 

(2T P - r t )(2(A - J) /T + 2r p - r t )) - 

8A tl A s2 (A/T + 2T P - T S )(J/T + T t - T s ) + 

16(A 2 tl A s2 + A 2 sl A t2 )(J/T + T t - T s ) 2 + 

2A t2 ({A/T + 2T p -J/T) 2 + 

r t (l + 2(J- A)/T-4T P )) + 

8A sl A t2 (T s ({A- J)/T-2T t ) + 

(J/T + T t )((J-A)/T + T t )- 

2T p (J/T + T t -T s ))+4A sl A tl x 

(8T 2 p + r a ((J - 2A)/T + r t ) - 2(A/T) x 

((j + a)/t + r t ) - 4r p ( j/t - t s + r t - 2A/t))\ 

for a strongly-coupled electron-phonon system in the 
bipolaronic regime, where we have introduced the dimcn- 
sionless parameters 



r« = 



^dEE 2 T a {E)df$/dE 
Tf~dEET a (E)dft/dE' 

J™dEE 3 T a (E)df$/dE 
T 2 j™dEET a (E)dfZ/dE' 



(43) 



(44) 
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and D\ t 2 = (1 



4.4,1,2 
Lp = 
L s ,t = 



-4A t i,2) _1 . Also 
(7p - r£), 

(is 



rL)/4. 



(45) 
(46) 



In the limit of a pure polaronic system (i.e. A t =A s 
0) the Lorenz number and Hall Lorenz number is 



L = L 



H 



Lp. 



If polarons are degenerate then L p = L 



(47) 
■ 2 /3 for 

any elastic scattering and any dimensionality. Indeed for 
a majority of metals L falls in the region of approximately 
(3.1 — 3.3) as expected in a degenerate Fermi liquid. How- 
ever Lp depends on the dimensionality and the scattering 
mechanism for non-degenerate polarons. In the opposite 
limit we have a pure singlet or triplet bipolaronic system 



(i.e. A t 



A r , 



or A s = A p = 0). We obtain from 



eqns.(41) and (42) 



L = Lh — L 



(48) 



respectively, which is about 6 times smaller than Lq. The 
last equation is expected for a charged Bose gas as noted 
by Mott and one of the authors (ASA)2&. In the general 
case our final equations, eqn.(41) and (42), yield Lorenz 
numbers that differ significantly from both limits at fi- 
nite temperatures. The main difference originates from 
the extra terms, which describe an interference of polaron 
and bipolaron contributions to the heat flow. In the low- 
temperature regime, T <C J and A, this contribution is 
exponentially small because the densities of triplet pairs 
and single polarons are small. However, this contribu- 
tion becomes important in the intermediate temperature 
range T c < T < T* because of the combination of the fac- 
tors (J/T) 2 and (A/T) 2 with the exponential form of the 
number densities. The contribution appears as the result 
of the recombination of a pair of polarons into triplet and 
singlet bound states at the cold end of the sample, which 
is reminiccnt to the contribution of the electron-hole pairs 
to the heat flow in semiconductors^. Our expressions 
for the Lorenz numbers eqns. (41,42) anticipate its values 
to deviate from Lq/4 for our system of predominantly 
charged Bose particles. Here the L terms are representa- 
tive of a type of scattering mechanism(s), which can be 
a combination of types. In the non-degenerate system it 
is given by 



(7- 



(r + 2) 



(49) 



where r is related to the energy dependence of the scat- 
tering time 



r cx E r . 



(50) 



VII. NUMERICAL ANALYSIS AND THE 
PHONON THERMAL CONDUCTIVITY 

Here we have shown that the present model fits the 
Hall Lorenz number, Lh, measured by Zhang et aln 




FIG. 3: Lh which fits nicely to experiment {YBa2Cus.Ofi.gs,) 
using Ob /Op = 0.36. This experiment showed a clear viola- 
tion of the W-F law. The inset gives the ratio of Hall Lorenz 
number to Lorenz number 



From these fits the phonon thermal conductivity can also 
be extracted. 

The charge and spin gap value was estimated by Mi- 
hailovic et ali2& for optimally doped Y BaiCu^O^x (x is 
the number of doped oxygen ions), giving A/2 = 675K 
and J = 150A", in their systematic analysis of charge and 
spin spectroscopies. According to Refi^i the main scat- 
tering mechanism above T c is the particle-particle colli- 
sions which gives a relaxation time T 8t t >p oc 1/T 2 . The 
chemical potential is pinned near the mobility edge, so 
that y — exp(/i/T) w 0.6 in a wide temperature range, 
if the number of localised states in the random poten- 
tial is about the same as the number of bipolarons^i. In 
Y BaiCuzO§+ x every excess oxygen ion x can localise a 
bipolaron so this approximation is reasonable. As a re- 
sult, there is only one fitting parameter in Lh, Eq.(42) 
which is the ratio of the bipolaron and polaron Hall 
angles 8;>/0p, where <d a = qBT a /m a and we assume 

e s = e t = e 6 . 

The model gives a good fit, as shown in Fig. 3, with a 
reasonable value of Ob/0 P = 0.36. By using the same 
single parameter for L, as was used to fit Lh to the ex- 
perimental data, we can see that the ratio of the Lorenz 
numbers varies with temperature with L being larger at 
lower temperatures and Lh growing at higher tempera- 
tures (inset in fig. 3). 

The model also describes the (quasi) linear in-plane 
resistivity, Hall resistivity T 2 -dependence and the Hall 
ratio, as observed in the cuprates (fig 4). As mentioned 
earlier the extraction of the phonon thermal conductiv- 
ity from experiment has proven difficult and inconclusive. 
Here we can extract this quantity by using our theoret- 
ically calculated L (lower inset of fig. 5) and using ex- 
perimental resistivity data™ (upper inset of fig. 5) . This 
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VIII. CONCLUSIONS 




300 



FIG. 4: T dependence of p xx , p xy and Rh for optimally doped 
cuprates. 



allows us to find K e i in the normal state, fig. 5. Using 
the same parameters as for the Lh fittings, L is seen to 
violate the W-F law (lower inset in fig. 5). One can see 
that the electronic thermal conductivity is therefore very 
weakly T-dependent. These results are similar to the 
findings by Takenaka et al and Salamon et al. To find 



0.012 




0.012 



0.011 



0.01 



0.009 



0.008 



0.007 



0.007 



300 



FIG. 5: Normal state K e l(= LT/p) in optimally doped 
cuprates calculated using the theoretical L (lower inset) and 
experimental resistivity data (upper inset). 

K p h we take for example Minami and Cooper— thermal 
conductivity data on YBa2Cuj,OQ.^ and subtract our 
K e i. The results are shown in fig. 6 implying that the 
T-dependence of k is predominately due to the phonon 
contribution. We believe that this method is the only way 
to extract n p h reliably. The lattice contribution to the 
diagonal heat flow appears to be much higher than it is 
anticipated in the framework of any Fermi-liquid model. 



Recent measurements by Proust et al^ 6 on 
Tl 2 Ba2Cu0e + s have suggested that the Wiedemann- 
Franz law holds perfectly well in the overdoped region 
and therefore conclude that the Fermi-liquid prevails at 
this doping. Alexandrov and Mott£, suggested that there 
might be a crossover from Bose-Einstein condensation to 
a BCS-like superconductivity across the phase diagram. 
Thus the results of Proust et al are compatible with 
the bipolaron picture. If the Fermi liquid does exist 
at overdoping then it is likely that the heavy doping 
causes an "overcrowding effect" where the polarons find 
it difficult to form bipolarons due to the larger number 
of competing holes 2 . Here we have shown that by the 



0.12 



0.11 



V 0.1 - 



0.09 



0.08 



0.07 



i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — 

Theoretical phonon thermal conductivity 

- -□ - - Experimental thermal conductivity 
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0.09 
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0.07 



100 
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200 
T(K) 



250 



300 



FIG. 6: Deduced and total k of Minami et al thermal 
conductivity data on Y Ba%CuzOe.&3 in the normal state. 

necessary inclusion of thermally excited polarons and 
triplets as the temperature rises, the bipolaronic gas can 
predict the Hall Lorenz number as found experimentally. 
The fits here have been on optimally doped cuprates, 
however it is most likely that this model will work in 
the underdoped as well. Also this analysis, along with 
experimental data, allowed us to estimate the phonon 
contribution to the thermal conductivity which has 
remained elusive in recent years. The interference of 
the polaron and bipolaron contributions to the energy 
flow breaks down the Wiedemann-Franz law and results 
in the unusual temperature dependence of the Lorenz 
numbers. This work further validates the Bipolaronic 
model and it is our belief that the superconducting state 
heat transport can be described by it as well. 

This work was supported by the Leverhulme Trust 
(grant F/00261/H) and by the EPSRC UK (grant 
R46977). We would like to thank Y. Zhang and PW. 
Anderson for helpful discussions of some experimental 
and theoretical aspects. 



9 



1 S. Chakravarty, R.B. Laughlin, D.K. Morr and C. Nayak, 
Phys. Rev. B, 63, 094503 (2001). 

2 A.S. Alexandrov and N.F. Mott, High Temperature Su- 
perconductors and Other Superfluids, (Taylor and Francis, 
London, 1994). 

3 V.J. Emery and S.A. Kivelson, Nature, 374, 434 (1995). 

4 P.W. Anderson, The Theory of Superconductivity in the 
High T c Cuprates, (Princeton University Press, Princeton, 
1997). 

5 S.C. Zhang, Science, 275, 4126 (1997). 

6 D. Pines, A.V. Balatsky and P. Monthoux, Phys. Rev. 
Lett, 67, 3448 (1991). 

7 G. Zhao, M.B. Hunt, H. Keller and K.A. Muller, Nature, 
385, 236 (1997). 

8 D. Mihailovic, CM. Foster, K. Voss and A.J. Heeger, Phys. 
Rev. B42, 7989 (1990). 

9 P. Calvani, M.Capizzi, S. Lupi, P. Maselli, A. Paolone, P. 
Roy, S-W Cheong, W. Sadowski and E. Walker, Solid State 
Commun. 91, 113 (1994). 

10 rp Timusk, C.C. Homes and W. Reichardt, in Anhar- 
monic properties of High T c Cuprates, eds. D. Mihailovic, 

G. Ruani, E. Kaldis and K.A. Muller, 171(World Scientific, 
Singapore, 1995). 

11 J.L. Cohn, S. Wolf and TA. Vanderah, Phys. Rev. B , 45, 
511 (1992). 

12 T. Egami, J. Low Temp. Phys. 105, 791(1996). 

13 A. Lanzara, P.V. Bogdanov, X.J. Zhou, S.A. Kellar, D.L. 
Feng, E.D. Lu, T. Yoshida, H. Eisaki, A. Fujimori,K. 
Kishio, J.I. Shimoyana, T. Noda, S. Uchida, Z. Hussain 
and Z.X. Shen, Nature, 412, 510 (2001). 

14 A. Chainani, T. Yokoya, T. Kiss, S. Shin, T. Nishio and 

H. Uwe, Phys. Rev. B, 64, 180509 (2001). 

15 G.M. Eliashberg, Sov.Phys. JETP, 11, 696 (1960). 

16 A.S. Alexandrov, Phys. Rev. B, 46, 2838 (1992). 

17 I.G. Lang and YA. Firsov, Sov. Phys. JETP, 16, 1301 
(1963). 

18 A.S. Alexandrov, Russ. J. Phys. Chem., 57, 273 (1983). 

19 A.S. Alexandrov and E.A. Mazur, Zh. Eksp. Teor. Fiz., 96, 
1773 (1989). 

20 J.G. Bednorz and K.A. Muller, Angew. Chem., Int. Ed. 
Engl. 57, 735 (1988). 

21 A.S. Alexandrov, Phys. Rev. B, 53, 2863 (1996). 

22 C.R.A. Catlow, M.S. Islam and X. Zhang, J. 
Phys.:Condensed Matter, 10 No. 3, L49 (1998). 

23 J. Bonca J and S.A. Trugman, Phys. Rev. B 64, 094507 
(2001). 

24 A.S. Alexandrov and P.E. Kornilovitch, J. 
Phys.:Condensed Matter, 14 No. 21, 5337 (2002). 

25 A.S. Alexandrov and N.F. Mott, J. Supercond (US), 7, 599 
(1994). 

26 J. Rossat-Mignod, L.P. Regnault, P. Bourges, C. Vettier, 
P. Burlet and J.Y. Henry, Physica Scripta, 45, 74 1992). 

27 H.A. Mook and M. Yethiraj, Phys. Rev. Lett, 70, 3490 
(1993). 

28 D. Mihailovic, V.V. Kabanov, K. Zagar, and J. Demsar, 



Phys. Rev. B60, 6995 (1999) and references therein. 

29 G. Deutscher, Nature, 397, 410 (1999). 

30 A.S. Alexandrov and A.F. Andreev, EuroPhys. Lett., 54, 
373 (2001). 

31 A.S. Alexandrov, A.M. Bratkovsky and N.F. Mott, Phys. 
Rev. Lett, 72, 1734 (1994). 

32 A.S.Alexandrov, A.M. Bratkovsky, N.F. Mott and E.K.H 
Salje, Physica C, 215, 359 (1993). 

33 K.A. Muller, G.M. Zhao, K. Conder and H. Keller, J. Phys. 
Cond. Mat., 10, L291 (1998). 

34 A.S. Alexandrov, Physica C, 305, 46 (1998). 

35 A.S. Alexandrov, Phys. Rev. B, 46, 14932 (1992). 

36 A.S. Alexandrov and N.F. Mott, Phys. Rev. Lett, 71, 1075 
(1993). 

37 A.S. Alexandrov, Doctoral thesis, Moscow Engineering 
Physics institute (1984). 

38 A.S.Alexandrov and C.J. Dent, Phys. Rev. B, 60, 15414 
(1999). 

39 A. Junod, Studies of High- Temperature Superconductors 
Vol. 19, eds. A. Narlikar, 1 (Nova science, Commack, NY, 
1996). 

40 A.S. Alexandrov, W.H. Beere, V.V. Kabanov and W.Y 
Liang, Phys. Rev. Lett., 79, 1551 (1997). 

41 A. Junod, A. Erb and C Renner, Physica C, 317-318, 333 
(1999). 

42 Y. Zhang, N.P. Ong, ZA. Xu, K. Krishana, R. Gagnon, 
and L. Taillefer, Phys. Rev. Lett., 84, 2219 (2000). 

43 K.K. Lee, A.S . Alexandrov and W.Y Liang, 
cond-mat/0211665 (2002). 

44 K. Takenaka, Y. Fukuzumi, K. Mizuhashi, S. Uchida, H. 
Asaoka and H. Takei, Phys. Rev. B, 56, 5654 (1997). 

45 R.C. Yu, M.B. Salamon, J. P. Lu and W.C. Lee, Phys. Rev. 
Lett., 69, 1431 (1992). 

46 R.W. Hill, C. Proust, L. Taillefer, P. Fournier and R.L. 
Greene, Nature, 414, 711 (2001). 

47 J. Takeya, Y. A ndo, S. Komiya and X.F. Sun, 
cond-mat/0108055 (2002). 

48 E. W. Carlson, V. J . Emery, S. A. Kivelson, and D. Orgad, 
cond-mat/0206217 and references therein. 

49 A. J. Leggett, Physica Fennica 8, 125 (1973); J Stat. Phys. 
93, 927 (1998); 

V. N. Popov, Functional Integrals and Collective Excita- 
tions (Cambridge: Cambridge University Press) (1987). 

50 A.S. Alexandrov, Physica C 363, 231 (2001). 

51 K.A. Muller, Phil. Mag. Lett, 82, 279 (2002). 

52 G.M. Zhao, Phys. Rev. B, 64, 024503 (2001). 

53 A. Anselm, Introduction of Semiconductor Theory , (Pren- 
tice and Hall, New Jersey, 1981). 

54 Kouji Segawa and Yoichi Ando, Phys. Rev. Lett., 86, 4907 
(2001). 

55 H. Minami and J.R. Cooper, Private communication. 

56 C. Proust, E. Boakin, R.W. Hill, L. Taillefer and A.P 
Mackenzie, |cond-mat/0202101| (2002) . 



